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Abstract 

The transition from the ordered commensurate phase to the incommen- 
surate gaussian phase of the antiferroelectric asymmetric six-vertex model is 
investigated by keeping the temperature constant below the roughening point 
and varying the external fields {h,v). In the {h,v) plane, the phase bound- 
ary is approached along straight lines 6v = k6h, where {Sh, 5v) measures the 
displacement from the phase boundary. It is found that the free energy singu- 
larity displays the exponent 3/2 typical of the Pokrovski-Talapov transition 
6f ~ const for any direction other than the tangential one. In the 
latter case 5f shows a discontinuity in the third derivative. 
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1 Introduction 



The asymmetric six-vertex model is the extension of the well-known symmetric 
six-vertex model to the case where external fields interact with the local fiuctuating 
variables ("dipoles" or "arrows") [|l], ||, |^. On a two-dimensional square lattice, two- 
valued variables live on links and interact at vertices, and, throughout this paper, 
we choose interactions that favor antiferroelectric order. An external homogeneous 
two-component field {h, v), favoring ferroelectric ordering, competes with the arrow- 
arrow couplings. 

If one treats the problem with the transfer matrix method, an exact solution 
is provided by the Bethe-ansatz ^, The phase diagram was outlined in the 
original paper [|l|, while further developments came more recently p|, |], |^, H (the 
list is not exhaustive. Papers dealing with the ferroelectric regime are not included). 
Since the early works appeared, it has been known that, in the {h,v) plane, there is 
a closed curve F separating an antiferroelectrically ordered (commensurate) massive 
region, inside F, from an incommensurate, disordered, massless region, outside F. 
The ordered region includes the symmetric model point (0, 0) and the free energy 
in it is field independent 

fh,u,h,v) = /o(7,m) 

The parameters 7, u contain the temperature dependence and will be defined in 
the next section. As F is approached from outside (i.e. from the massless phase), 
keeping 7 and u fixed, the free energy is expected to display a singularity with a 
characteristic exponent 3/2. More precisely, the following result has been proven in 

i- 

Let {h{b) ,v{b)) be a parametric equation of F, to be given later, with b some real 
parameter running in a bounded interval. Lieb and Wu found that, if /i = 0, which 
implies 6 = 0, the singular part of the free energy is 

/(7, u, 0, U(0) + 6v) = /o(7, u) + 6f{j, u, 6v) = /o(7, u) + 0(7, u) {Svf/' (1) 

yielding a divergence ~ 5v^^^'^ for the susceptibility. (A divergence with exponent 
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1/2 was also found for the specific heat at fixed field) . Borrowing a terminology in- 
troduced in later years, the transition belongs to the Pokrovski-Talapov universality 
class [|, 

Eq.(|l|) does not clarify how the singularity depends on both field components. 
A further step in this direction was recently achieved by Noh and Kim ^ who 
extended the method of ||rT| to relate macroscopic quantities such as susceptivities 
to finite size corrections. They showed that, in the critical phase 



9V 

dh'^ dhdv 
dvdh dv"^ 
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where g is the coupling constant of the gaussian model on which the critical incom- 
mensurate phase renormalizes. As F is approached from this phase, g —>■ 2. 

In this paper, the question is settled of how, at fixed 7, u the singularity of f{h, v) 
depends on both {h,v). While finding the leading singularity as a function of {h,v) 
simultaneously has proven to be elusive, it will be shown that the singularity does 
depend on the direction in the {h, v) plane. 

The results are best summarized if one introduces the following notation. Set 

A = vthi- htvi 

Next, consider a variation into the incommensurate phase h = h{b)+6h, v = v{b)+6v 
where 6v = kSh and k fixes a slope not tangential to F. Then 

6f{6h) = {^y/^^[2{k-^-i)6hf^ (3) 
A Stt ht 

so the singularity is governed by an exponent 3/2 for all these directions. Yet, if F 
is approached tangentially, which amounts to take 5v = j^Sh 

, c+6h^ if 6h>0 
6f{6h) = { (4) 

c.5h^ if 5h<0 



where c+ 7^ c_ are 6-dependent and are given in Eq. (pS]). In other words there is a 
jump in the third derivative. The calculation breaks down an the two points on F 
where ht = 0, i.e. where the tangent to F is parallel to the f-axis. Those cases were 
examined in |^ and an analogous conclusion was reached. 

The paper is divided as follows. Section 2 presents a summary of the Bethe- 
ansatz and already known results. Section 3 deals with a perturbative expansion 
of the integral equations typical of the Bethe-ansatz and section 4 exploits that 
expansion to determine the singularities of the free energy. 



2 Basic definitions and summary of Bethe-ansatz 

The Boltzmann weights of the six allowed configurations are grouped into RZ{n) 
as shown in Fig.l. In the framework of the Bethe-ansatz, the spectral parameter 
notation, chosen here, seems more natural than the traditional one. Row-to-row 
transfer matrices 

N 

/3 k=l 

with periodic boundary conditions {Pn+i = Pi) commute for different values of the 
spectral parameter 

[T{u),T{u')]=0 

Arrow conservation at each vertex, and periodic boundary conditions, imply that 
T{u) breaks into blocks between states with the same number of up (and down) 
arrows. Let n be the number of arrows reversed with respect to the reference state 
I If, . . . , |>. The Bethe-ansatz provides the following solution to the eigenvalue 
problem for T{u) i, i |, i: 

^,(^_2„)+hivrSinh(7 -u)^^{\ sinh(| + m - 



A(U) = e^iN-2n)+hNr ^^^^^^KI ^J .N JT 



sinh7 ^-J-^ smh{^ - u + '-^) 
^,(iv_2n)-fejVr SinhM ^jv sinh(-f + m - ^) 



sinh7 f-Ji sinh(f -u + ^ 
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is an eigenvalue if the "rapidities" {aj}, j = 1, 2, . . . , n satisfy the set of equations 



Sinh(f + _r M„2hN A ^^Ml + " 



Wh(i-^)^ ^ ' Msinh(7-tK-a.)) ' ''"'•••'^ 

In the hmit oo, the rapidities {a^} condense into curves in the complex plane, 
and are conveniently described by a density function R{q) 

27r 

Riaj) = lim — 

^ jv^oo N{aj+i - aj) 

Eqs. (^) are replaced by a single linear integral equation that governs the thermody- 
namics of the model. Introduce the functions (this notation is somewhat redundant, 
but it has been adopted in many previous papers and it will be kept here) 

,sinh(f + f ) sinh(7+f). 



p'{a) = -t\og[ . 'I I' e(a) = -zlog[- 



sinh(f-f) ' ' ^^sinh(7-f)- 



da da 
and the vertical polarization 

2v 

y = hm (1 - -) (7) 

N—i-oo iV 

Then, for a state described by a rapidity curve C, the density /2(a;) solves |^, 

[ df3K{a - P)R{P) = R{a) (8) 

and 

p%a) - — / ci/3e(a - + 2ih = 2nx - ^—^ <x< ^—^ (9) 

where x is the real parameter of the curve. Let A = —a + ib, B = a + ib he the 
two endpoints of the curve. (We take for granted that B = —A* because we wish 
to consider, in each sector of fixed n, the largest transfer matrix eigenvalue, which 
is real and unique by Perron- Froboenius theorem ||13|. Since {aj} —>■ {—a*} is a 
symmetry of (M) we expect it to hold for the solution corresponding to the unique 
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largest eigenvalue in each sector). Solution of (|^) implicitly depends on the endpoints 
A, B and contributes to make h dependent on a, h through 0, |, |] 

l-y 1 fB 



2 

1 



— [ daR{a;A,B) (10) 
2tt J a 



p\A)+p\B)- — I dl3R{l3-A,B)[Q{A-p) + Q{B~l3)]+Aih = {) (11) 
In J A 

In the transfer matrix formalism, the free energy is determined by the largest eigen- 
value Aq, so, neglecting a factor 

/(m, 7, h,v) = - hm (12) 

N—*oo iV 

With an abuse of language, the relevant eigenstate will be called "ground state" . 
From (I), all eigenvalues are such that A(m) = Aji{u) + Al{u). One of the two 
addends dominates over the others for specific values of the parameters. We set 

lim ^ In Aij(M) = Fii{u, 7, h, y) + vy = 

= h + ln ^'""^^^ '""^ l- daR(a; A, B)fR(a, u) (13) 
smh 7 2tt J a 

lim ^lnAL(M) = FL{u,'y,h,y)+vy = 

siiili u \ 

= -h + \n — \ / daR(a:A,B)fL(a;u) (14) 

smh 7 2iT J A 



with 



sinh(i + ^-f) sinh(-f + «-f) 

fR{a;u) = \n——^ fL{a;u) = \n- 



sinh(f-u + f) ' ^ sinh(^-M + f) 

and call F = max{FR,Fi}. The equilibrium polarization and the free energy are 
then determined by 

f{u,'y,h,v)= mm {-F{u,-f,h,y)-vy)= min J^{u,-f,h,v,y) (15) 
— i<j/<i — i<?/<i 

The ground state solution of (|) is explicitly computed when h and v are sufficiently 
small |l|, ^ 1^, It corresponds to n = N/2 and, in the N —>■ 00 limit, to a = tt, 
—7 < 6 < 7. Then (|^) can be solved by Fourier transform. Even though the 
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solution has appeared many times in the hterature, it is worthwhile to recall it here 
to introduce the elliptic function notation 

R{a- -vr + ^6, vr + ^h) = g ""J = :^dn(:^; k) (16) 

n±^oo2cosh(n7) vr vr 

where = ^ and I{k)(I'(k)) is the complete elliptic integral of the first kind with 
modulus k{k' = Vl - k^). Replacing (0) into (|10|) and (0) yields 0, |, |, g 

y = h{b) = Z{-h) (17) 

where 

z{x) ^ + V tr^j^pi (18) 

2 ^ n cosh(n7) 
Within the range —7 < 6 < 7 there is a crossing between A/j and A^,. Let's call 

id the point where the ground state curve meets the imaginary axis in the a-plane. 

It turns out that A^j dominates when d < 7 — 2m while A^ does when d > 7 — 2u. 

One can check it by replacing ([161) into (|T3| ) , (|T^ ) , or, more quickly, by the following 

argument. Select chains such that n = y is an odd number. By virtue of the 

symmetry {aj} — » {—a*} one a, say a, has to be pure imaginary. The other n — 1 

{aj} can be paired to give a positive contribution to A^ and A^, so the contribution 

of a determines the sign of A.r and A^. It is easily seen that 

AR(a) > (<)0 if Ima < (>)7 - 2u 
AL(a) > (<)0 if Ima > (<)7 - 2u 

In the thermodynamic limit Ima d, and since by Perron-Froboenius theorem 
AAT^max niust bc positivc, we get the desired result. When crossing, though, A^j and 
Al connect smoothly, so no singularity of / appears and one finds from (p!3D,(p^ 
the field independent value of the symmetric six-vertex model |]T^ 

+00 g~27n 

f{u, 7, h, t>) = — 2 — — - sinh(nu) sinh 72(7 — u) < u < 7 (19) 

n cosh(7n) 

The region in the (/i, v) plane where ([T9|) is valid is bounded by (0) and the value 
of the w-field at which the ground state moves away from the y = Q sector. This is 



fixed by the equation 

dF 

V{h) = --^Ih fixed, y=0 (20) 

that yields i |J, i 

U(6) = Z(7- |7-2n-6|) -7<6<7 (21) 

More precisely, v{b) in (^) runs over half of the curve, the other half is obtained 
from the symmetry f{h,v) = f{—h, —v), which follows trivially from reversal of all 
arrows in the statistical sum. Eqs. (|T7|) and (|21| ) provide the parametric equation 

(h{b),v{b)) of r. 

3 The expansion in 6a, 6b 



Eqs. (|]), (|TOD, (|TTD , ([T3|), ([T^ ) determine the thermodynamics of the system. In 



principle, one should solve (^ and plug R{a; A, B) into the others. This would give 

y = y{a,b) h = h{a,b) F = F{a,b) 

Once the first two are inverted and plugged into the third, one can minimize 
— F(7, u, h, y) — vy w.r. to y keeping h, v fixed. Of course, (|) is not exactly 
solvable for a ^ it. Therefore we attempt an expansion that generalizes the method 
that Lieb and Wu applied to the h = case, where the rapidities are real (6 = 0) 
and everything is a function of a only. Namely, it will be assumed that a unique 
solution exists for (|^) at least in a narrow neighborhood of the segment a = vr, 
-7 < 6 < 7 and that the dependence of R{a;A,B), y{A,B), h{A,B), F{A,B) is 
analytical on A and B in this neighborhood. This assumption is partly warranted 
by the fact that, if a solution exists for (^, it develops a pole at a = ±Z7, where 
C,{<y) has poles. In the following calculation, the endpoints are kept far away from 
the singularities of the inhomogeneous term. With this assuption, variations can be 
computed by taking derivatives in the relevant integral equations. One starts with 
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dARia;A,B) + ^ ['^ dl3K{a- (3)dARi(3;A,B) = ^ K {a - A) R{A; A, B) (22) 

dBR{a;A,B) + — dpK{a- p)dBRiP;A,B) = — K (a - B)R{B; A, B) (23) 
2tx J a 2% 

and so forth by taking further derivatives. When A = Aq = — vr + ib and B = Bq = 
n + ib, —7 < 6 < 7 all these equations can be solved by Fourier transform. For 
instance 

„ +00 nb-j\n\ 

dAR{a;Ao,Bo) = 77^ E e-*"°(-)"- = -dsRoia; A^, B^) 

2vr„f^oo 2 cosh 717 

ci = R{Ao;Ao,Bo) = R{Bo;Ao,Bo) 

Notice that ci = if 6 = ±7. This case was dealt with in p|. The solution is 
replaced into the analogous expansion for y{A,B), h{A,B) and Ffi{A,B) (we will 
consider only the domain where > F^, that is d < 7 — 2u). As an example, at 
the first order 

dv 2 2 

51 = -2ii 3^R(o,AB)da + -R(A;A,B) 

dv 2 2 

JL = dBR{a;A,B)da-—R{B;A,B) 
oB 2tt J a 27C 

For the field, after using (|) and ( ]TT| ) 

dh 1 1 

-4^^ = df3dARW;A,B)[e{A-P) + eiB-P)] + R{A-A,B)[l + —e{B-A)] 

OA Zn J A ZTT 

dh 1 1 

-4^7^^ = -TT/ dpdBR{P;A,B)[Q{A~P) + Q{B~P)] + R{B-A,B)[l + —Q{B-A)] 

and finally 

^ = dadAR{a;A,B)fn{a;u)-—R{A;A,B)fn{A-u) + — 

oA 2tt J A 271 oA 

dF 1 1 9/1 

a# = i^j^dadBR{a-A,B)fR{a-u) + —R{B-A,B)fR{B-u) + — 

To evaluate the last two equations it is useful to know the branch cut structure of 
u). The cuts run from 2(7 — 2u) to +ioo and from —iipf + 2-u) to — icx). The 
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transition from the "A/j regime" to the "A^ regime" occurs when the integration 
path crosses the branch cut at i{'y — 2u), since u) has a branch cut also starting 
from i{'~f — 2u) but running all the way down to —ioo. Furthermore Im//j(Ao; u) = n 
and lmfji{Bo; u) = -tx. 

It is now obvious how to go on by taking derivatives. The expansion has been 
carried out up to the third order. The variation of each quantity is a third order 
polynomial in 5a, 6b. It is convenient to introduce a more compact notation that 
brings out the geometrical meaning of the coefficients involved. Define 

hJb) = —h(b) = }{-) ; = ^dn(/ H : k) 

^ ' db ^ ' 2 1 coshn7 vr ^ ti ' ' 



where we have used ([16|) to express the series at hand as elliptic functions. For v{b) 
one has to distinguish between the two cases 6>7 — 2nor6<7 — 2u, but the 
elliptic function expression is the same for both cases 

v^ib) = ^v{b) = ^dn(/ + -i2u + b)-k) 
ab IT IT 

We will also need 

, /,N dht .,o,Hk).o ,^ ilb ,^ ilb 
ab IT 7T vr 

v,{b) = ^ = -tk\-^fsn{I+-{b + 2u);k)cn{I+-{b + 2u);k) 

db 71 TV 71 

ht is negative in the interval —7 < 6 < 7, and vanishes when b = ±7. At these two 
points, the tangent to F is parallel to the f-axis. Instead, vt = when b = 7 — 2m. At 
this point the tangent to F is parallel to the h-axis. There is of course another point 
of this kind on the other half of F, obtained by the symmetry {h,v) {—h, —v). 
Finally it is proven in appendix A that the combination vthi —viht, that will appear 
later, is definite negative. For polarization and horizontal field one has 

6y = ^6a - ^6a' + ^6a6b + - ^)5a' - ^5aHb + ^5a5b' (24) 

^ TT 27r TT 27r ^ 2 3 ^ 27r 27r ^ ' 

6h = h,6b - ^6a^ - ^6a6b - ^6a^ - ^6a^6b + ^6b' (25) 
2 2 6 2 6 ^ ^ 
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whereas, if we set = Z{'-f — I7 — 2m — b\)/27T 

6F = —2htZo5a + htX^ZoSa^ — {2hiZo H — —)5a6b 

+ ^ h Zo(-3/itXo + 2x3) — + (— h hiXoZojda 6b 

- + ^ + Zox,)5a6b' (26) 

Here and X3 are non-zero terms that play httle role in what follows 
^0 = -(I + 2EA— ) 

TT ^ cosh 77-7 

= -E(-) — T — = i^):r^M^^-.k)U=±n+ib 

^ coshn7 vr da"^ vr 

What actually has to be minimized though is JF = —F — vy. Using (|2T|) and ( ^6]) 

5jF = — — v{b)5y — 5v5y = —5v6y H — —5a6b 

71 

- ( + )^ + ( + + — — 6a 6b (27) 

TT TT O TT 27r 2tT 

These expansions reduce to those of @| when b ±7 (the limit should be taken in 
the elliptic functions because several series are not convergent when |6| = 7). No 
term 56" appears in (0), and that had to be expected since the line a = vr (i.e. 
6a = 0), —7 < 6 < 7 corresponds to y = 0. Furthermore, while (§) and following 
equations make perfectly sense for any values of A and B, one might object that 
solutions of d^) can always be taken to have — vr < Rea < vr. Hence a should vary 
in [0,7r]. The question at issue is what solutions (|^) admit when n > N/2, knowing 
that, when n = N/2, they already fill a line stretching from — vr to vr. Fortunately, 
the question can be bypassed. Only variations 6v > will be considered and it is 
physically clear that 6v > tends to align arrows "up", therefore brings 6y > 0, 
that is n < N/2 and, from (P^, 6a < 0. This is sufficient because only the upper 
half of r, given by (|l3), (l2l| ) is being considered. The other half, where to drive the 
system into the incommensurate phase one needs a variation 6v < 0, can, as usual, 
be recovered from f{h,v) = f{—h, —v). 
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4 Minimization of T and free energy singularity 



The minimum of T should now be taken with respect to when 7, m, /i, v are kept 
fixed. Keeping bh fixed at a given value means that ba and bh are not independent. 
Two difi^erent ways will be followed to deal with (|2^ ) , (|25|) and (|27|) and they give the 
same results. 

Neglecting terms 5a^, b}? in (EBI) and solving for bh one finds 



55 = ^ - ^5a2 - ^bh^ + 7^5a5/i + 0(5a^ bc?bh, babh^ , dh") (28) 
/Z^ 2:1X1 ^'^t 2rii 

when inserted into one gets 



5^ = ba{-^bv-l^bvbh+'^bh+'^^^^l±^bh^ + ...) 



+ 5a2(^^5^ + C25/.2 + ...) + ^A + ... (29) 

/vr DVT 

A = - (30) 

All terms neglected are higer order, i.e. 5a^, ba^bh, ba?bvbh'^ , babvbh"^, etc. and it will 
soon become clear that dropping them is justified. The coefficient C2 is 6-dependent 
and to know its specific value will not be necessary in the following. 

Let's consider the coefficient of ba in (^). It is easy to see, from ([T7| ) and 
(pT|) that it vanishes when the variation {bh, dv) is taken along F. This had to 
be expected. The curve F is described by minima of T falling on the line a = tt, 
—7 < & < 7 or, stated otherwise, with ba = 0. In fact, the vanishing of the first 



order term in (29) implies that bJ-" has a stationary point at ba = 0. Suppose next 



that we approach F along any direction other than the tangential one, that is 

bv = kbh ky^^ (31) 

ht 

To move into the incommensurate phase one has to take bh > for k > vt/ht and 
bh < for k < Vt/hf. Clearly the linear terms in the coefficient of ba in ( p9D do not 
cancel and therefore the terms bvbh, bh"^ (and higher) in the same coefficient can be 
dropped. Same thing can be said about the ba"^ term since, obviously, ba"^ << ba. 
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The term 5a^ must be retained in all cases, because its coefficient is finite. Instead 
terms like 5a?5h, 5a^ etc. are clearly negligible compared with the 5a^ term. It can 
be seen by inspection that the terms thrown away in solving (|25| ) are also negligible. 
The whole procedure assumes though that a small variation 5v) brings about a 
small displacement ((5a, 5h) in the minimum of T. The upshot is that it is legitimate 
to keep 

6T 



— A + (5a( — -5v 

DTT TT 



TT 



Sh) 



(32) 



Once (|31| ) is taken into account, the solution of 



occurs at 



6ac] 



d5a 

'2{kht-vt)5h 
A 



1/2 



(33) 



Notice that the previous limits on k guarantee that the solution is real. The sign 
has been chosen to make sure that 5a < 0. An elementary check shows that 



d^{5T)/d{5af 
into 



> and it confirms that 5ao is a minimum. When 5ao is plugged 
one finds that, when 5v = k6h, k ^ 

6f = /(7, u,h + Sh,v + 5v) - /o(7, u) = [hfnh.[2{k - ^)5hf' (34) 

So the exponent is 3/2 for all these directions. We pass next to the case where the 
transition line is approached tangentially, i.e. 

ht 



We have then 



= A + ^^6h6a^ 

Gtt 27r 



6a6h^ 



2Tlhl 



The minimum lies at 

(5ao = - 



ht 



/A 



(35) 



(36) 
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Since d"^ {5J-') / d{5aY\saQ > 0, 6ao is indeed a minimum. Inserting it into 6J^ yields a 
jump in the third derivative of the free energy. Namely, at htSv = vtSh 

, c+6h^ if 6h>0 
6f = f{^,u,h + 6h,v + 6v)-fo{i,u) = { ^ (37) 



c.5h^ if Sh<0 



where 



1 



Ci = -—^{9±J9^ + A^){g' + 2A'±gJg^ + A^) g = VtXtX, (3 



When 2m + 6 = 7, f( = 0. Then ( P5| ) somewhat simplifies, since g = Q and the 
variation tangential to V is just 5h, so 

An analogous result, with 5v ^ 6h was reached in but only for the two points 
b = ±7, where the tangent to F is parallel to the f-axis 

A second method of calculation will now be presented. From (|15]), the equation 
determining the equilibrium value of y is 

dF 

^U + . = (39) 

Eq. (pOD is just a particular case of this condition. It identifies the value of v for 
which the minimum occurs at y = 0. If one thinks to invert (|15|), (|2^), F becomes 
a function of a, b, so 

-.-^1^1 ,^1^1 ft- ft am 

oa uy uo uy 

which makes sense if ( plf ) and (|25|) are indeed invertible, that is if the denominator 
in the RHS does not vanish. It is immediate to see from (p^), (p5|) that this is true 
for the points we are considering. Instead, the points b = ^'y{ht = 0) are saddle 
points for h{a, b) p|. Expanding the denominator in (^0]) and retaining terms up to 



the second order ( to have the third order in (^0]) one should expand y, h and F to 
the fourth order) 

—V = —2'kZq — Vt6b H — -60^ -Sb"^ '^SaSb 

" 2 2 2 
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Writing v = v{b) + 6v, owing to (^) one arrives at 

Sv = vt6b-^{6a^ -6b^) + ^6a6b + ... (41) 
Sh = ht6b--{6a^-6b^)-^6a6b + ... (42) 

whose solution yields 6a{6h, Sv) and Sb{Sh, 6v). Their result has to be used in (|^) to 
produce the leading singular part of 6f. Again, two cases have to be distinguished. 
If 6v = k6h, k 7^ Vt/ht, consistency of (|40|), ( ^T]) requires that 6b ~ c6a^ + 0(60^), so 
at the leading order 

Sv = {cvt - 

5h = {cht — :r)6a^ 



Hence 



1 t>i — khx 

c = — — oa = — 

2 Vt — kht 



2{kht-vt)6h 



A 

which coincides with (pS]). 6 J-" is obviously the same previously used, so the free 
energy coincides with (^). Suppose instead 6v = k6h but k = Vt/hf. Clearly (0), 
(0) give 6b = 6h/ht at the leading order and 

6h'^ 

Ma^ + 2vtXo6a6h - A— 5- = 

hf 

whose solution is, again (PB|). Yet, replacing 6b = 6h/ht and ( PB| ) into ( P7[ ) is not 
correct. The reason is that a term 6a6b appears in ( pTf ) that would require solving 
(PI) and ( ^21 ) for 56 up to the next order in 6h. The problem can be bypassed by 



replacing e.g. (^) into the term 6a6b of (^). So one gets 

= - — A + ^^^^^6h6a? — ^^^—^6h6a6b + {hiVt + htVi) 
Gvr 27r tt/i^ 27r 

and after setting 56 = 6h/ht one is back to ( P^D and, because of (PU]), to the 5/ of 
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A Appendix 



One has to prove that A < for all b in (—7,7). Setting 



x = I + 



ilb 



TT 



from (H) 



A = {-YikMn{x; k)dn{y; k) 

TT 



y = I+—{h + 2u) 

TT 



sn{x] k)cn{x; k) sn(?/; k)cn{y; k) 
dn(x; k) dn(?/; k) 



Perform Landen's trasform Ill4 
k ki - 
sn(Mi;/ci) : 



1-k' 

1 + k' 
k sn(u; k)cn{u] k) 



M — > Ml = (1 + k')u 



kl'^ dn{u; k) 



The new half-periods are 



i[ = i'{k,) = (1 + k')i'{k) 



and so 



A = kl^hk{-f dn{x; k)dn{y; k) [sn(^(6 + 2n); ki) - sn(^; ki)] 
vr 7 7 



A table of signs of dn, cn and sn is given in ||TJ]. By inspecting the possible cases 
one sees that A < always. A possible exception comes from dn{y\ k) = that is 
when 6 + 2u = 7. In this case Vt = and 

A = -viht{b = 7 - 2u) = htk'{-f < 



TT 



There might be a simpler proof. 
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^h+v sinh(7-M) _ ^f3(s/2-e)+h+v 



sinh7 



2 — e 



-h-v sinh(7-M) _ ^p(^ij2-e)-h-v 
sinh7 



2 — e 



-h+v sinhu 
sinh7 



= e 



l3(-S/2-e)-h+v 



r^h—v sinhM 
sinh 7 



^Pi-5/2-e)+h-v ^ ^ll(^) 



Rll{u) 



Figure 1: Boltzmann weights in the notation with spectral parameter u compared 
to that of ref . [2] . The physical region is < m < 7 
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